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The use of weak measurements for performing quantum tomography is enjoying increased atten¬ 
tion due to several recent proposals. The advertised merits of using weak measurements in this 
context are varied, but are generally represented by novelty, increased efficacy, and foundational 
significance. We critically evaluate two proposals that make such claims and find that weak mea¬ 
surements are not an essential ingredient for most of their advertised features. 


I. INTRODUCTION 

The business of quantum state tomography is convert¬ 
ing multiple copies of an unknown quantum state into 
an estimate of that state by performing measurements 
on the copies. The naive approach to the problem in¬ 
volves measuring different observables (represented by 
Hermitian operators) on each copy of the state and con¬ 
structing the estimate as a function of the measurement 
outcomes (corresponding to different eigenvalues of the 
observables). Though tomography can be performed in 
such a way, there are more general ways of interrogat¬ 
ing the ensemble; indeed, generalizations such as ancilla- 
coupled [1] and joint [2] measurements lead one to eval¬ 
uate the problem of tomography from the perspective 
of generalized measurements [3], an approach which has 
yielded many optimal tomographic strategies [4-8]. 

An interesting subclass of generalized measurements is 
the class of weak measurements [9-13]. Figure 1 gives 
a quantum-circuit description of a weak measurement. 
Weak measurements are often the only means by which 
an experimentalist can probe her system, thus making 
them of practical interest [14-20]. Sequential weak mea¬ 
surements are also useful for describing continuous mea¬ 
surements [21]. 

Weak measurements are also central in the more con¬ 
tentious formalism of weak values [22]. In particular, the 
technique of weak-value amplifieation [23] has generated 
much debate over its metrological utility [24-34]. 

The two proposals we review in this paper assert that 
it is useful to approach the problem of tomography with 
weak measurements holding a prominent place in one’s 
thinking. Some care needs to be taken in identifying 
whether a particular emphasis has the potential to be 
useful when thinking about tomography, given the large 


body of work already devoted to the subject. Since weak 
measurements are included in the framework of general¬ 
ized measurements, none of the known results for optimal 
measurements in particular scenarios are going to be af¬ 
fected by shifting our focus to weak measurements. In 
Sec. II we outline criteria for evaluating this shift of fo¬ 
cus. 



FIG. 1: A circuit depicting an ancilla-coupled measurement. 
Here A is a two-system Hermitian operator, j?/;) is the state 
of the system being measured, |0) is the initial state of the 
meter, e is a real number parameterizing the strength of the 
measurement, and O is a standard observable with outcomes 
Oj. If |e| <C I the measurement is weak, U{e) ^ 1, and very 
little is learned or disturbed about the system by measuring 
the meter. 


We apply these criteria to two specific tomographic 
schemes that advocate the use of weak measurements. 
Direet state tomography (Sec. Ill) utilizes a procedure 
of weak measurement and postselection, motivated by 
weak-value protocols, in an attempt to give an oper¬ 
ational interpretation to wavefunction amplitudes [35]. 
Weak-measurement tomography (Sec. IV) seeks to out¬ 
perform so-called “standard” tomography by exploiting 
the small system disturbance caused by weak measure¬ 
ments to recycle the system for further measurement [36] . 


II. EVALUATION PRINCIPLES 
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tomography. The primary tool we utilize is generalized 
measurement theory, specifically, describing a measure¬ 
ment by a positive-operator-valued measure (POVM). A 
POVM assigns a positive operator Ep to every measur¬ 
able subset F of the set ft of measurement outcomes 
X G 12. For countable sets of outcomes, this means the 
measurement is described by the countable set of positive 
operators, 

{^x}xen • (^•^) 

The positive operators Ep are then given by the sums 

Ef=EEx- ( 2 - 2 ) 

xeF 

For continuous sets of outcomes the positive operator as¬ 
sociated with a particular measurable subset F is given 
by the integral 

Ep= [ dE^ . (2.3) 

Jf 

These positive operators capture all the statistical prop¬ 
erties of a given measurement, in that the probability of 
obtaining a measurement result y within a measurable 
subset F C O for a particular state p is given by the 
formula 

Pr(x e F\p) = Tr {pEp) ■ (2.4) 

That each measurement yields some result is equivalent 
to the completeness condition, 

En = t. (2.5) 

POVMs are ideal representations of tomographic mea¬ 
surements because they contain all the information rele¬ 
vant for tomography, i.e., measurement statistics, while 
removing many irrelevant implementation details. If 
two wildly different measurement protocols reduce to 
the same POVM, their tomographic performances are 
identical. 


A. Novelty 

The authors of both schemes we evaluate make claims 
about the novelty of their approach. These claims seem 
difficult to substantiate, since no tomographic protocol 
within the framework of quantum theory falls outside 
the well-studied set of tomographic protocols employing 
generalized measurements. To avoid trivially dismissing 
claims in this way, however, we define a relatively conser¬ 
vative subset of measurements that might be considered 
“basic” and ask if the proposed schemes fall outside of 
this category. 

The subset of measurements we choose is composed of 
randomly chosen one-dimensional orthogonal projective 
measurements [hereafter referred to as random ODOPs; 


see Fig. 2(a)]. These are the measurements that can be 
performed using only traditional von Neumann measure¬ 
ments, given that the experimenter is allowed to choose 
randomly the observable he wants to measure. This is 
quite a restriction on the full set of measurements allowed 
by quantum mechanics. Many interesting measurements, 
such as symmetric informationally complete POVMs, like 
the tetrahedron measurement shown in Fig. 2(b), cannot 
be realized in such a way. With ODOPs assumed as 
basic, however, if the POVM generated by a particular 
weak-measurement scheme is a random ODOP, we con¬ 
clude that weak measurements should not be thought of 
as an essential ingredient for the scheme. 

Identifying other subsets of POVMs as “basic” might 
yield other interesting lines of inquiry. For example, 
when doing tomography on ensembles of atoms, weak 
collective measurements might be compared with non- 
adaptive separable projective measurements [15, 20]. 

B. Efficacy 

Users of tomographic schemes are arguably less inter¬ 
ested in the novelty of a particular approach than they 
are in its performance. There is a variety of performance 
metrics available for state estimates, some of which have 
operational interpretations relevant for particular appli¬ 
cations. Given that we have no particular application 
in mind, we adopt a reasonable figure of merit, Haar- 
invariant average fidelity, which fortuitously is the figure 
of merit already used to analyze the scheme we consider 
in Sec. IV. This is the fidelity, /(p, p{x)) ? of estimated 
state p{x) with the true state p, averaged over possible 
measurement records y and further averaged over the 
unitarily invariant (maximally uninformed) prior distri¬ 
bution over pure true states. For the case of discrete 
measurement outcomes, this quantity is written as 

F{p,E) := f dpY,P^ix\p)fip,Kx)) ■ (2.6) 

X 

An obvious problem with this figure of merit is its de¬ 
pendence on the estimator p. We want to compare mea¬ 
surement schemes directly, not measurement-estimator 
pairs. To remove this dependence we should calculate 
the average fidelity with the optimal estimator for each 
measurement, expressed as 

F{E) := maxF(p,E) . (2.7) 

P 

To avoid straw-man arguments, it is also important 
to compare the performance of a particular tomographic 
protocol to the optimal protocol, or at least the best 
known protocol. Both proposals we review in this pa¬ 
per are nonadaptive measurements on each copy of the 
system individually. Since there are practical reasons for 
restricting to this class of measurements, we compare to 
the optimal protocol subject to this constraint. 
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FIG. 2: (a) Implementation of a random ODOP by performing a randomly selected [probability P(A)], basis-changing unitary 
U{X) before making a projective measurement of a standard observable O, with outcomes oj. (b) POVM elements, represented 
as Bloch vectors, for the tetrahedron measurement, an example of a POVM that cannot be implemented as a random ODOP, 
because the POVM elements cannot be sorted into sets of equally weighted orthogonal projectors. 


This brings up an interesting point that can be 
made before looking at any of the details of the weak- 
measnrement proposals. For our chosen figure of merit, 
the optimal individual nonadaptive measurement is a 
random ODOP (specifically the Haar-invariant measure¬ 
ment, which samples a measurement basis from a uniform 
distribution of bases according to the Haar measure). 
Therefore, weak-measurement schemes cannot hope to 
do better than random ODOPs, and even if they are 
able to attain optimal performance, weak measurements 
are clearly not an essential ingredient for attaining that 
performance. 


C. Foundational significance 

Many proposals for weak-measurement tomography 
are motivated not by efficacy, but rather by a desire to 
address some foundational aspect of quantum mechan¬ 
ics. This desire offers an explanation for the attention 
these proposals receive in spite of the disappointing per¬ 
formance we find when they are compared to random 
ODOPs. 

There are two prominent claims of foundational sig¬ 
nificance. The first is that a measurement provides 
an operational interpretation of wavefunction amplitudes 
more satisfying than traditional interpretations. This is 
the motivation behind the direct state tomography of 
Sec. Ill, where the measurement allegedly yields expec¬ 
tation values directly proportional to wavefunction am¬ 
plitudes rather than their squares. 

The second claim is that weak measurement finds a 
clever way to get around the uncertainty-disturbance 
relations in quantum mechanics. The intuition behind 
using weak measurements in this pursuit is that, since 


weak measurements minimally disturb the system being 
probed, they might leave the system available for further 
use; the information obtained from a subsequent mea¬ 
surement, together with the information acquired from 
the preceding weak measurements, might be more infor¬ 
mation in total than can be obtained with traditional 
approaches. Of course, generalized measurement theory 
sets limits on the amount of information that can be ex¬ 
tracted from a system, suggesting that such a founda¬ 
tional claim is unfounded. We more fully evaluate this 
claim in Sec. IV. 


III. DIRECT STATE TOMOGRAPHY 

In [35] and [37] Lundeen et al propose a measurement 
technique designed to provide an operational interpre¬ 
tation of wavefunction amplitudes. They make various 
claims about the measurement, including its ability to 
make “the real and imaginary components of the wave- 
function appear directly” on their measurement device, 
the absence of a requirement for global reconstruction 
since “states can be determined locally, point by point,” 
and the potential to “characterize quantum states in situ 
... without disturbing the process in which they feature.” 
The protocol is thus often characterized as direct state 
tomography (DST). 

To evaluate these claims, we apply the principles dis¬ 
cussed in Sec. IF Lundeen et al. have outlined procedures 
for both pure and mixed states. We focus on the pure- 
state problem for simplicity, although much of what we 
identify is directly applicable to mixed-state DST. To 
construct the POVM, we need to describe the measure¬ 
ment in detail. The original proposal for DST of Lundeen 
et al. calls for a continuous meter for performing the weak 
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measurements. As shown by Maccone and Rusconi [38], 
the continuous meter can be replaced by a qubit meter 
prepared in the positive cr^ eigenstate |+), a replacement 
we adopt to simplify the analysis. Since wavefunction 
amplitudes are basis-dependent quantities, it is necessary 
to specify the basis in which we want to reconstruct the 
wavefunction. We call this the reconstruction basis and 
denote it by {|^)}o<n<d’ where d is the dimension of the 
system we are reconstructing. 

The meter is coupled to the system via one of a collec¬ 
tion of interaction unitaries where 

(3.1) 

The strength of the interaction is parametrized by (p. A 
weak interaction, i.e., one for which \p\ 1, followed by 

measuring either (jy or cTz on the meter, effects a weak 
measurement of the system. In addition, after the in¬ 
teraction, there is a strong (projective) measurement di¬ 
rectly on the system in the conjugate basis 
which is defined by 


{n\cj) = LO^gVd , CO := . (3.2) 

The protocol for DST of Lundeen et al., motivated by 
thinking in terms of weak values, discards all the data 
except for the case when the outcome of the projective 
measurement is cq. This protocol is depicted as a quan¬ 
tum circuit in Fig. 3. 


1+) 

w 



(<Ty, <7z^ - i 
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FIG. 3: Quantum circuit depicting direct state tomogra¬ 
phy. The meter is coupled to the system via one of a family 
of unitaries, {F(^,n}o<n<d? each of which corresponds to a 
reconstruction-basis element. The meter is then measured in 
either the y or z basis to obtain information about either the 
real or imaginary part of the wavefunction amplitude of the 
selected reconstruction-basis element. This procedure is post- 
selected on obtaining the co outcome from the measurement 
of the system in the conjugate basis. While the postselection 
is often described as producing an effect on the meter, the 
circuit makes clear that the measurements can be performed 
in either order, so it is equally valid to say the measurement 
of the meter produces an effect on the system. 


measurement, are given by 

2 sinp 


(c 


vnn,co dPr(co|C„,^,^) 


Re(V’nT*) 


sm2p — 2 sinp 2 


(3.3) 


d PT{co\Un,cp,1p) 


I'lpn 


{^^)\n,ca ^ 

where 'ijn •= 
come Co is 


2 sinp 




(3.4) 


Pr (cq \ 7/^) 

The probability for obtaining out 


Pr(co|t/n,c^,V') 

= + 2(cos</?- l)[Re('*/)„T*) - IV’nP]) , 


(3.5) 


and 


T;=^V’n- (3.6) 

n 

We can always choose the unobservable global phase of 
lip) to make T real and positive. With this choice, which 
we adhere to going forward, {o'y)\^ provides informa¬ 
tion about the real part of -0^, and {cyz)\nco Pi*ovides 
information about the imaginary part of 

Specializing these results to weak measurements gives 

= ^(K)lco,n + *('^^)lco,n) +C>(<^^). (3.7) 

This is a remarkably simple formula for estimating the 
state \ 2 p)\ There is, however, an important detail that 
should temper our enthusiasm. Contrary to the claim 
in [37], this formula does not allow one to reconstruct the 
wavefunction point-by-point (amplitude-by-amplitude in 
this case of a finite-dimensional system), because one has 
no idea of the value of the “normalization constant” T 
until all the wavefunction amplitudes have been mea¬ 
sured. This means that while ratios of wavefunction 
amplitudes can be reconstructed point-by-point, recon¬ 
structing the amplitudes themselves requires a global re¬ 
construction. Admittedly, this reconstruction is simpler 
than commonly used linear-inversion techniques, but it 
comes at the price of an inherent bias in the estimator, 
arising from the weak-measurement approximation, as 
was discussed in [38]. 

The scheme as it currently stands relies heavily on 
postselection, a procedure that often discards relevant 
data. To determine what information is being discarded 
and whether it is useful, we consider the measurement 
statistics of ay and az conditioned on an arbitrary out¬ 
come Cm of the strong measurement. To do that, we 
first introduce a unitary operator Z, diagonal in the re¬ 
construction basis, which cyclically permutes conjugate- 
basis elements and puts phases on reconstruction-basis 
elements: 


For each n, the expectation values of ay and cr^, condi¬ 
tioned on obtaining the outcome cq from the projective 


^ Wj) — ki+i) 5 


Z\n) =uj^\n) . (3.8) 
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As is illustrated in Fig. 4, postselecting on outcome Cm 
with input state \ 2 p) is equivalent to postselecting on cq 
with input state l^p) = 1^)* 



(b) 

FIG. 4: (a) Postselection on outcome Cm, achieved by post- 
selecting on Co after cyclic permutation of the conjugate basis 
by application of the unitary Z~^. (b) Since Z commutes 
with (a) is identical to postselecting on co with input 

state Z~^ |'0). 

Armed with this realization, we can write reconstruc¬ 
tion formulae for all postselection outcomes, 

■ ( 3 - 9 ) 

This makes it obvious that all the measurement outcomes 
in the conjugate basis give “direct” readings of the wave- 
function in the weak-measurement limit. Postselection 
in this case is not only harmful to the performance of 
the estimator, it is not even necessary for the interpre- 
tational claims of DST. Henceforth, we drop the post¬ 
selection and include all the data produced by the strong 
measurement. 

The uselessness of postselection is not a byproduct 
of the use of a qubit meter. In the continuous-meter 
case, the conditional expectation values in the weak- 
measurement limit are given as weak values 

Weak-value-motivated DST postselects on meter out¬ 
come p = 0 to hold the amplitude {p\x) constant and thus 
make the expectation value proportional to the wave- 
function {xl^p). Since {p\x) is only a phase, however, it is 
again obvious that postselecting on any value of p gives a 
“direct” reconstruction of a rephased wavefunction. Shi 
et al. [42] have developed a variation on Lundeen’s pro¬ 
tocol that requires measuring weak values of only one 


meter observable. This is made possible by keeping data 
that is discarded in the original postselection process. 

We now consider whether the weak measurements in 
DST contribute anything new to tomography. It is al¬ 
ready clear from Eqs. (3.3) and (3.4) that for this pro¬ 
tocol to provide data that is proportional to amplitudes 
in the reconstruction basis, the weakness of the interac¬ 
tion is only important for the measurement of ay. We 
are after something deeper than this, however, and to 
get at it, we change perspective on the protocol of Fig. 3, 
asking not how postselection on the result of the strong 
measurement affects the measurement of ay or az, but 
rather how those measurements change the description 
of the strong measurement. As is discussed in Fig. 3, 
this puts the protocol on a footing that resembles that of 
the random ODOPs in Fig. 2(a). 

The measurement of az , which provides the imaginary- 
part information, is trivial to analyze, because the analy¬ 
sis can be reduced to drawing more circuits. In Fig. 5(a), 
the interaction unitary is written in terms of system uni- 
taries Un,± •= that are controlled in the 2 ;- 

basis of the qubit. The az measurement on the meter 
commutes with the interaction unitary, so using the prin¬ 
ciple of deferred measurement, we can move this mea¬ 
surement through the controls, which become classical 
controls that use the results of the measurement. The 
resulting circuit, depicted in Fig. 5(b), shows that the 
imaginary part of each wavefunction amplitude can be 
measured by adding a phase to that amplitude, with the 
sign of the phase shift determined by a coin flip. This 
is a particular example of the random ODOP described 
by Fig. 2(a). We conclude that weak measurements are 
not an essential ingredient for determining the imaginary 
parts of the wavefunction amplitudes. 

Measuring the real parts is more interesting, since the 
ay measurement does not commute with the interaction 
unitary. We proceed by finding the Kraus operators that 
describe the post-measurement state of the system. The 
strong, projective measurement in the conjugate basis 
has Kraus operators Km = \cm){cm\^ whereas the unitary 
interaction followed by the measurement of ay with 
outcome ±, has (Hermitian) Kraus operator 

:= {±y\U^,n 1+) 

= T(^i + (y25^_i)l„)(n|), (3.11) 

s± := ± sm{p ± 7r/4), 


where the eigenstates of ay are |±^) := |0) + 

e^*7^/4 |1) )/\/2. The composite Kraus operators, 
yield POVM elements := 

\cm){cm\K^±"^. For each n, 
these POVM elements make up a rank-one POVM with 
2d outcomes. 
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FIG. 5: (a) The imaginary-part measurement in a circuit 
where the interaction unitary is written in terms of system 
unitaries Un,± •= controlled in the z-basis. This 

circuit is equivalent to that of Fig. 3 when cr^ is measured on 
the meter; the result of the measurement reveals the sign of 
interaction, i.e., which of the unitaries Un,± was applied to the 
system, (b) The coherent controls can be turned into classical 
controls by moving the measurement before the controls. The 
result is a particular instance of the random ODOP depicted 
in Fig. 2(a). 


The POVM elements can productively be written as 


Tp{y,n) 

^±,m 

(y) 

= 

Ay,n)\/ Ay,n) 
^±,m lY±,m 

5 

(3.12) 

dy,'n)\ 

^±,m / 

V 

II 

|c„) 



— 

) + (p\/2s± — . 

L)a;™” 1 

n)) / \l 2 Ap , 





(3.13) 

(y) 



(3.14) 

POVM for each n does not 

fit into 

our framework of 


random ODOPs, but can be thought of as within a wider 
framework of random POVMs. Indeed, the Neumark 
extension [39, 40] teaches us how to turn any rank-one 
POVM into an ODOP in a higher-dimensional Hilbert 
space, where the dimension matches the number of out¬ 
comes of the rank-one POVM. 

Vallone and Dequal [41] have proposed an augmen¬ 
tation of the original DST to obtain a “direct” wave- 
function measurement without the need for the weak- 
measurement approximation. The essence of their pro¬ 
tocol is to perform an additional ax measurement on 
the meter. The statistics of this measurement allow the 
second-order term in cp to be eliminated from the real- 
part calculation, giving a reconstruction formula that is 


exact for all values of (p. Of course, the claim that the 
original DST protocol “directly” measures the wavefunc- 
tion is misleading, and directness claims for Vallone and 
DequaFs modifications are necessarily more misleading. 
Even ratios of real parts of wavefunction amplitudes no 
longer can be obtained by ratios of simple expectation 
values, since these calculations rely on both ax and ay 
expectation values for different measurement settings. 

We analyze this additional meter measurement in the 
same way we analyzed the ay measurement. The Kraus 
operators corresponding to the meter measurements are 

^ ('cos(/:? - l) \n){n\ , 

(3.15) 

j^(x,n) _ ^_| ^ smip\n){n\ . (3.16) 

The composite Kraus operators, = KmK±'^\ 

yield POVM elements that can be written as 



It is useful to ponder the form of the POVM elements 
for the ay and ax measurements of the DST protocols. 
For the original DST protocol of Fig. 3, without post¬ 
selection, the only equatorial measurement on the me¬ 
ter is of ay] the corresponding POVM elements, given 
by Eq. (3.12), are nearly measurements in the conjugate 
basis, except that the n-component of the conjugate ba¬ 
sis vector is changed in magnitude by an amount that 
depends on the result of the ay measurement. Eor the 
augmented DST protocol of Vallone and Dequal, the ad¬ 
ditional POVM elements (3.17), which come from the 
measurement of ax on the meter, are quite different de¬ 
pending on the result of the ax measurement. Eor the 
result +, the POVM element is similar to the POVM el¬ 
ements for the measurement of cr^, but with a different 
modification of the n-component of the conjugate vector. 
Eor the result —, the POVM element is simply a measure¬ 
ment in the reconstruction basis; as we see below, the ad¬ 
dition of the measurement in the reconstruction basis has 
a profound effect on the performance of the augmented 
DST protocol outside the region of weak measurements, 
an effect unanticipated by the weak-value motivation. 

Although claims regarding the efficacy of DST are 
rather nebulous, we consider the negative impact of the 
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weak-measurement limit on tomographic performance. 
In doing so, we assume for simplicity that the system 
is a qubit, in some unknown pure state that is specified 
by polar and azimuthal Bloch-sphere angles, 0 and (j). 
In this case we assume that the reconstruction basis is 
the eigenbasis of the conjugate basis is the eigenbasis 
of ax. 

The method we use to evaluate the effect of variations 
in ip is taken from the work of de Burgh et al. [43], 
which uses the Cramer-Rao bound (CRB) to establish 
an asymptotic (in number of copies) form of the average 
fidelity. 



FIG. 6: (Color online) CRB C{p) of Eq. (3.22) for origi¬ 
nal DST (solid black) and augmented DST (dotted green for 
probability 1/4 for ax and ay measurements and 1/2 proba¬ 
bility for a az measurement; dashed red for equal probabilities 
for all three measurements). As the plot makes clear, the op¬ 
timal values for p are far from the weak-measurement limit. 
Values of p for which ~ 0 give exceptionally large CRBs, 
confirming the intuition that weak measurements learn about 
the true state very slowly. The CRB for original DST also 
grows without bound as p approaches 7r/2, since that mea¬ 
surement strength leads to degenerate Kraus operators and 
a POVM that, not informationally complete, consists only of 
projectors onto ax and ay eigenstates of the system. The 
CRB remains finite when the meter measurements are aug¬ 
mented with (7x, since the resultant POVM oX p — tt!2 then 
includes az projectors on the system [i.e., projectors in the re¬ 
constructions basis; see Eq. (3.19)], giving an informationally 
complete overall POVM. 

In analyzing the two DST protocols, original and aug¬ 
mented, we assume that the two values of n are chosen 
randomly with probability 1/2. For the original protocol, 
we choose the az and ay measurements with probability 
1/2. For the augmented protocol, we make one of two 
choices: equal probabilities for the ax^ ay, and az mea¬ 
surements or probabilities of 1 /2 for the az measurement 
and 1/4 for the ax and ay measurements. Formally, these 
assumed probabilities scale the POVM elements when all 
of them are combined into a single overall POVM. 


The asymptotic form involves the Fisher informations, 
Je and J^, for the two Bloch-sphere state parameters, 
calculated from the statistics of whatever measurement 
we are making on the qubit. The CRB already assumes 
the use of an optimal estimator. When the number of 
copies, N, is large, the average fidelity takes the simple 
form 

FP) ^ 1 - , (3.21) 

C{p) = f dO sinO 
Jo 

(3.22) 

1/ 1 sin^ 6> \ 

^ Jo 

Though we have derived analytic expressions for the 
Fisher informations, it is more illuminating to plot 
the CRB C(p), obtained by numerical integration (see 
Fig. 6). For original DST, the optimal value of p is just 
beyond 7r/4, invalidating all qualities of “directness” that 
come from assuming p 1. For the augmented DST of 
Vallone and Dequal, the optimal value of p moves toward 
7r/2, even further outside the region of weak measure¬ 
ments. In both cases, C{p) blows np at p = 0; for weak 
measurements, C{p) is so large that the information gain 
is glacial. 

We visualize this asymptotic behavior by estimating 
the average fidelity over pure states as a function of N 
using the sequential Monte Carlo technique [46] , for var¬ 
ious protocols and values of p. Figure 7 plots these re¬ 
sults and shows how the average fidelity, for the optimal 
value of p, approaches the asymptotic form (3.21) as N 
increases. We note that the estimator used in these sim¬ 
ulations is the estimator optimized for average fidelity 
discussed in [6] . If we were to use the reconstruction for¬ 
mula proposed by Lundeen et al, the performance would 
be worse. 

Our conclusions are the following. First, postselection 
contributes nothing to DST. Its use comes from attention 
to weak values, but post selection is actually a negative 
for tomography because it discards data that are just 
as cogent as the data that are retained in the weak-value 
scenario. Second, weak measurements in this context add 
very little to a tomographic framework based on random 
ODOPs. Finally, the “direct” in DST is a misnomer [47] 
because the protocol does not provide point-by-point re¬ 
construction of the wavefunction. 

The inability to provide point-by-point reconstruction 
is a symptom of a general difficulty. Any procedure, 
classical or quantum, for detecting a complex amplitude 
when only absolute squares of amplitudes are measurable 
involves interference between two amplitudes, say, A and 
B, so that some observed quantity involves a product 
of two amplitudes, say, Re(A*5). If one regards A as 
“known” and chooses it to be real, then Re(5) can be 
said to be observed directly. This is the way amplitudes 
and phases of classical fields are determined using inter¬ 
ferometry and square-law detectors. 













measurements on the system. The optimality claim can¬ 
not be true, of course, since a random ODOP based on 
the Haar invariant measure for selecting the ODOP basis 
is optimal when average fidelity is the figure of merit, but 
the novelty of the information extraction remains to be 
evaluated. 



Qi 
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FIG. 7: (Color online) Average infidelity 1 — F as a function of 
the number N of system copies for three measurements. The 
dashed red curve is for augmented DST with equal probabil¬ 
ities for the three meter measurements; the value ip — 1.25 is 
close to the optimal value from Fig. 6. The other two curves 
are for original DST: the solid black curve is for ip = 0.89, 
which is close to optimal (this curve nearly coincides with 
the dashed red curve for augmented DST); the dashed-dotted 
purple curve is for a small value p = 0.1, where the weak- 
measurement approximation is reasonable. The three dotted 
curves give the corresponding asymptotic behavior C{p)/N. 
The two weak-measurement curves illustrate the glacial in¬ 
formation acquisition when weak measurements are used; the 
dashed-dotted curve hasn’t begun to approach the dotted 
asymptotic form for N = 10^. 


Of course, quantum amplitudes are not classical fields. 
One loses the ability to say that one amplitude is known 
and objective, with the other to be determined relative to 
the known amplitude. Indeed, if one starts from the to¬ 
mographic premise that nothing is known and everything 
is to be estimated from measurement statistics, then A 
cannot be regarded as “known.” DST fits into this de¬ 
scription, with the sum of amplitudes, T of Eq. (3.6), 
made real by convention, playing the role of A. The 
achievement of DST is that this single quantity is the 
only “known” quantity needed to construct all the am¬ 
plitudes 'ipn from measurement statistics. Single quantity 
or not, however, T must be determined from the entire 
tomographic procedure before any of the amplitudes 
can be estimated. 


IV. WEAK-MEASUREMENT TOMOGRAPHY 

The second scheme we consider is a proposal for qubit 
tomography by Das and Arvind [36]. This protocol was 
advertised as opening up “new ways of extracting infor¬ 
mation from quantum ensembles” and outperforming, in 
terms of fidelity, tomography performed using projective 


FIG. 8: Quantum circuit depicting the weak-measurement to¬ 
mography protocol of Das and Arvind. Two identical meters 
are used as ancillas to perform the weak and x measure¬ 
ments. The circuit makes clear that there is nothing impor¬ 
tant in the order the measurements are performed after the 
interactions have taken place, so we consider the protocol as 
a single ancilla-coupled measurement. 

The weak measurements in this proposal are measure¬ 
ments of Pauli components of the qubit. These measure¬ 
ments are performed by coupling the qubit system via an 
interaction unitary, 

[/(i) = g-i<T,0F^ (4.1) 

to a continuous meter, which has position Q and momen¬ 
tum P and is prepared in the Gaussian state 

The position of the meter is measured to complete the 
weak measurement. The weakness of the measurement is 
parametrized by e = 1/Ag^. 

The Das-Arvind protocol involves weakly measuring 
the 2 ) and x Pauli components and then performing a 
projective measurement of ay. We depict this protocol 
as a circuit in Fig. 8. Das and Arvind view this protocol 
as providing more information than is available from the 
projective ay measurement because the weak measure¬ 
ments extract a little extra information about the 2 : and 
X Pauli components without appreciably disturbing the 
state of the system before it is slammed by the projective 
ay measurement. Again, we turn the tables on this point 
of view, with its notion of a little information flowing 
out to the two meters, to a perspective akin to that of 
the random ODOP of Fig. 2(a). We ask how the weak 
measurements modify the description of the final pro¬ 
jective measurement. For this purpose, we again need 
Kraus operators to calculate the POVM of the overall 
measurement. 
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FIG. 9: (Color online) Effective distributions over measurement bases visualized as distributions over the positive-y Bloch 
hemisphere. Very weak measurements of z and x (e.g., e = 0.1) don’t perturb the final y measurement very much, so the 
distribution of bases is concentrated about the y axis. Very strong measurements (e.g., e = 1) end up extracting most of the 
information in the first z measurement, so the distribution becomes concentrated around the z axis. The optimal measurement 
(e ^ 0.575) has an effective distribution that is nearly uniform. 


The Kraus operators for the projective measurement vectors, 
are = |±^)(±^|, and the (Hermitian) Kraus op¬ 
erator for a weak measurement with outcome q on the G{qi^q 2 ) 

meter is 



e(gi + gj + 2) ^ 


K^i\q) := {q\ |(/.) 


X cosh eqi cosh eq 2 , 


(4.6) 





(4.3) 


X cosh(eg/2) + cTj sinh(eg/2)^ 

The Kraus operators for the whole measurement pro¬ 
cedure are ^^±(^ 1 ,^ 2 ) •= K^^\q 2 )K^^\qi). From 

these come the infinitesimal POVM elements for out¬ 
comes gi, g 2 , and ±: 

dE±{qi,q2) 

■■= Kl{qi,q2)K±{qi,q2) (4.4) 

= K^^\q,)K^^\q 2 ) \±y){±y\K^^\q 2 )K^^\qi) . 

These POVM elements are clearly rank-one. 

Using the Pauli algebra, we can bring the POVM ele¬ 
ments into the explicit form, 

dE±{qi,q2) = i^l(gi, g2)-f^±(9i, 92 ) 

= dqi dq2 G{qx,q2)\{t + n±(gi, 92) • ct) , 

( 4 . 5 ) 

where we have introduced a probability density and unit 


n±(9i,g2) 


X sinh eq2 ± y + z sinh egi cosh eq2 
cosh eqi cosh eq2 


(4.7) 


We note that G(gi,g 2 ) = G(—gi,—^ 2 ) and n±(gi,g 2 ) = 
—fi=p(—gi, —g 2 ). This means that the overall POVM is 
made up of a convex combination of equally weighted 
pairs of orthogonal projectors and is therefore a random 
ODOP. From this perspective, the weak measurements 
are a mechanism for generating a particular distribution 
from which different projective measurements are sam¬ 
pled; i.e., they are a particular way of generating a dis¬ 
tribution P(A) in Fig. 2. Several of these distributions 
are plotted in Fig. 9. 

It is interesting to note that the value of e that Das and 
Arvind identified as optimal (about 0.575) produces a 
distribution that is nearly uniform over the Bloch sphere. 
This matches our intuition when thinking of the mea¬ 
surement as a random ODOP, since the optimal random 
ODOP samples from the uniform distribution 

To visualize the performance of this protocol, we again 
use sequential Monte Carlo simulations of the average 
fidelity. Das and Arvind compare their protocol to a 
measurement of da,, cr^, and cr^, whose eigenstates are 
mutually unbiased bases (MUB). In Fig. 10, we compare 
Das and Arvind’s protocol for e = 0.575 to a MUB mea¬ 
surement and to the optimal projective-measurement- 
based tomography scheme, i.e., the Haar-uniform ran¬ 
dom ODOP. 
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FIG. 10: (Color online) Average infidelity 1 — F as a function 
of the number N of system copies for three measurements: 
Das and Arvind’s measurement protocol (dotted-dashed blue) 
with e = 0.575; MUB consisting of Pauli ax, ay^ and az mea¬ 
surements (solid black); random ODOP consisting of projec¬ 
tive measurements sampled from the Haar-uniform distribu¬ 
tion (dashed red). The dotted lines, 1/N and 13/12iV, are 
the CRBs dehned by Eq. (3.22) for the optimal generalized 
tomographic protocol and MUB measurements, respectively. 

We don’t discuss the process of binning the position- 
measurement results that Das and Arvind engage in, 
since such a process produces a rank-2 POVM that is 
equivalent to sampling from a discrete distribution over 
projective measurements and then adding noise, a prac¬ 
tice that necessarily degrades tomographic performance. 

We conclude that the protocol does not offer anything 
beyond that offered by random ODOPs and that its claim 
of extracting information about the system without dis¬ 
turbance is not supported by our analysis. In particular, 
when operated optimally, it is essentially the same as the 
strong projective measurements of a Haar-uniform ran¬ 
dom ODOP. It is true that the presence of the 2 : and 
X measurements provides more information than a pro¬ 
jective y measurement by itself, but this is not because 
the 2 : and x measurements extract information without 
disturbing the system. 

V. SUMMARY AND CONCLUSION 

Our analysis of weak-measurement tomographic 
schemes gives us guidance for future forays into tomog¬ 


raphy. 

POVMs contain the necessary and sufficient informa¬ 
tion for comparing the performance of tomographic tech¬ 
niques. Specific realizations of a POVM might provide 
pleasing narratives, but these narratives are irrelevant for 
calculating figures of merit. Optimal POVMs for many 
figures of merit and technical limitations are known. A 
new tomographic proposal should identify restrictions 
on the set of available POVMs that come about from 
practical considerations and compare itself to the best 
known POVM in that set. The question of the opti¬ 
mality of Das and Arvind’s tomographic scheme is easily 
answered by identifying what POVMs arise from “projec¬ 
tive measurement-based tomography” and realizing these 
POVMs are optimal even in the generalized nonadaptive, 
individual-measurement scenario. 

Claims about novel properties of a state-reconstruction 
technique should be evaluated as a comparison with a 
motivated restriction on the set of POVMs. The false di¬ 
chotomy between “tomographic methods” and whatever 
new method is being proposed obfuscates that all new 
methods implement a POVM and that reconstructing a 
state from POVM statistics is nothing but tomography. 
Our analysis shows that even the relatively bland and 
conceptually simple set of random ODOPs captures most 
of the behavior exhibited by more exotic protocols. 

It is appropriate to move beyond the minimal, 
platform-independent POVM description when consid¬ 
ering ease of implementation or when trying to provide 
a helpful conceptual framework. Nonetheless, a pleasing 
conceptual framework should not be confused with an 
optimal experimental arrangement. If the experimental 
setup described by Lundeen et al. happens to be the eas¬ 
iest to implement in one’s lab, the state should still be 
reconstructed using techniques developed in the general 
POVM setting rather than the perturbative reconstruc¬ 
tion formula presented in work on DST, regardless of how 
attractive one finds the wavefunction-amplitude analogy. 
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